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Abstract. In a recent publication, J. Phys.: Condens. Matter 14 13777 (2002), 
Kuzovkov et al announced an analytical solution of the two-dimensional Anderson 
localisation problem via the calculation of a generalised Lyapunov exponent using 
signal theory. Surprisingly, for certain energies and small disorder strength they 
observed delocalised states. We study the transmission properties of the same model 
using well-known transfer matrix methods. Our results disagree with the findings 
obtained using signal theory. We point to the possible origin of this discrepancy and 
comment on the general strategy to use a generalised Lyapunov exponent for studying 
Anderson localisation. 
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It is generally believed that in the absence of both spin-orbit scattering and 
magnetic fields, there is no metal-insulator transition for non-interacting electrons 
in disordered two-dimensional (2D) systems. This belief is based on the scaling 
hypothesis [T] and supported by detailed finite-size scaling analysis of numerical data 

Nevertheless, an analytical proof that there are no extended states in a disordered 
2D system described by the Anderson model is still missing. 

Recently, Kuzovkov et al studied disordered 2D systems using signal theory, 
and they found parameter regimes where the system under consideration transmits the 
incident signal. They interpret this result as an indication for a metallic phase. This 
strongly contradicts standard wisdom, and, therefore, in this Comment, we reconsider 
the transmission properties of a 2D Anderson model using well-known transfer matrix 
methods. Similar calculations were presented before by Pendry 

Consider a 2D Anderson model on the M x L lattice (M is the width and L is the 
length of the system). The discrete Schrodinger equation 

may be rewritten using the transfer matrix T^^^ as 

Cr)^ ..'■>(:•: J, ^ 

Here, = Ho + en is the (M x M) Hamiltonian of the nth slice which contains random 
(uncorrelated) energies and \E'„ is the vector (■i/'ni, '^n2, • • • , ipnu)- For simplicity we 
assume for the disorder potentials that (£„m) = and {Snm^n'm') = (^"^^nn'^mm' ■ The 
angle brackets denote averages over disorder. 

In order to study the properties of quantities {'ipnm^'nm) (signals) as considered by 
Kuzovkov et al [3^ we construct the tensor product of Eq. (0), 

(:r)«(:r)-^^'(tJ«(:L)- 

The size of the matrix T^^) = j-a) (g, ^a) ig AM"^ x AM"^. We may now average Eq. © 
over the disorder. Since the matrix T^^^ and the tensor product on which it operates 
are statistically independent, we obtain a non-random matrix T^'^^ = (T^^-'), 

/ (T^l (g) 1 + Do ® Do -Do ® 1 -l®Do 1 \ 



rp{2) 



(4) 



Do 1 (D -10 

1(g) Do -1 

V 1 / 

with Dq = E — Hq. The matrix Q can be transformed by the matrix Q Q, where Q 
diagonalises Hq. Using Dirichlet boundary conditions in the transversal direction, one 
obtains 

TT 

QDqQ'^ = K, Ki = E-2coski, ki = _^ -^ h i = l,...,M. (5) 



Comment on Kuzovkov et al 



3 



As a result, the 4M^ eigenvalues A of T^^^ may be calculated from 



M 



Hdet 



/ KiKj + Cr^ — A —Ki —Kj 1 \ 
Ki -A -1 
-1 -A 



0. (6) 



V 1 -A ; 

Therefore, the eigenvalues fulfil the equations 

- + or^) + + 1^] - 2) - /(/ti/tj- - a^) + 1 = 0. (7) 

Using / = exp(i(/), the Eqs. ((Zj) read 

2 cos 2q — 2/€jKj cos + (k^ + — 2) = 2cr^z sin q. (8) 

For the ordered case (cr = 0) the 4M^ solutions are obtained as Qij = ±aj ± with 
2 costtj = E — 2 cos fcj, and the eigenvalues of the unperturbed matrix are given by 

l{a^ = 0) = exp [i{±ai ± aj)] . (9) 

If cr 7^ 0, then Eq. (jSJ shows that q is real only (i.e. |A| = 1), if g is zero or a multiple of 
TV. Therefore, the only eigenvalues of the transfer matrix T^'^^ which lie on the unit circle 
are / = ±1. Note that these eigenvalues are independent of the disorder. Therefore, 
it is evident that they correspond to some internal symmetry of the Anderson model 
(see also Ref. In fact, from Eqs. ((Tj) one finds M eigenvalues A = 1 (one for each 
i = j). They just correspond to current conservation in each of the M channels j3]. 
We will call these eigenvalues 'trivial eigenvalues'. Eigenvalues / = —1 are obtained for 
corresponding to the special energy E = cos /cj + cos kj, which can be discarded 
as irrelevant. 

Our results so far do not disagree with those of Ref. 0. However, Kuzovkov et 
al claim that in the limit M oo eigenvalues |A| = 1 will appear, which do not 
correspond to the trivial eigenvalues discussed above. Mathematically, the limit M — oo 
of the discrete model discussed here bears various conceptual and technical difficulties. 
Therefore we studied this limit numerically and find no new solutions with |A| = 1 apart 
from the trivial ones. It therefore appears that the metallic solutions found by Kuzovkov 
et al are an artefact of a questionable mathematical limit procedure. 

In conclusion, our analysis shows that apart from the 'trivial' eigenvalues / = ±1, 
which correspond to symmetries of the 2D Anderson Hamiltonian, all other eigenvalues 
have an absolute value different from unity. This holds for any M. Therefore, all input 
signals are damped away if the system is sufficiently long. This result contradicts the 
findings in Ref. 0. We believe that the reasons for this contradiction are as follows: 
Kuzovkov et al reduced the size of the transfer matrix from 4M^ x 4M^ to only ~ M. 
This drastic reduction was achieved by introducing a translational symmetry in the 
transversal direction. A precise description and motivation of this additional averaging 
procedure is absent in their paper. We do not believe that such a symmetry really 
holds in general. We therefore have to conclude that by this reduction of the size of 
the transfer matrix they transformed the model defined in Eq. into another model. 
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But even then Kuzovkov et al do not find eigenvalues |A| = 1 different from tlie trivial 
eigenvalues. Only using a questionable limit procedure M ^ oo such eigenvalues are 
obtained. 

Finally, we would like to comment on the general strategy to tackle the Anderson 
localization problem using methods similar to the one discussed here and in Ref . |S] . We 
do not think that the analysis of the second moments of the wave function l-^p (or of any 
higher moments as was done here as well as in Ref. is able to detect a metallic 

phase. The physical reason is the following: In order to detect a metallic phase one 
needs to average the logarithm of the wave function itself [H]. The procedures discussed 
here, while mathematically correct, always lead to quantities which are damped away 
for large systems, i.e. yield eigenvalues |A| 7^ 1 apart from the trivial ones. This can be 
exemplified most drastically with an analysis of a 3D system: Here, one obtains Eq. (jH)) 
with Ki = E — 2 cos ka — 2 cos kp, i = M{a — 1) + f3, a, P = 1, . . . , M. Consequently, 
the absolute values of all non-trivial eigenvalues differ from unity independently of the 
strength of the disorder. Since weakly disordered 3D systems are metallic, we conclude 
that in contrast to (ln(|?/^|)) an analysis of (l^/'P) does not help to observe a metallic 
regime even if it exists. 
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